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CN ■ Abstract 

^ ■ We propose an analytical study of relativistic tunneling through opaque barriers. We obtain a 

closed formula for the phase time. This formula is in excellent agreement with the numerical 
simulations and corrects the standard formula obtained by the stationary phase method. An 
important result is found when the upper limit of the incoming energy distribution coincides with 
the upper limit of the tunneling zone. In this case, the phase time is proportional to the barrier 
width. 

I.INTRODUCTION 

. The question of how long it takes a particle to tunnel through a barrier potential has been the issue of 

tIh ' intriguing discussions in the last decades. Since the first studies regarding the tunneling process [1,2], 

several tunneling time definitions have been proposed [3-5]. However, there is no general agreement 
about a satisfactory definition. Considering a time-dependent description in terms of wave packets, 
^ ■ Hartman [6] used the stationary phase method (SPM), previously applied to scattering problems [7,8], 

to estimate the instant in which the transmitted peak appears in the free region after the barrier. This 
time definition is known as phase time. One can also introduce a time integration over the probability 
^ of finding the particle inside the barrier. This average time spent in the potential region, regardless 

whether transmission or reflection occurs, is known as dwell time [9, 10]. Phase and dwell times have 
a well established mutual relation [3, 11]. The important point to be noted here is that, in the opaque 
limit, both predict the independence of tunneling times on the barrier width, the so-called Hartman 
effect [6] . It is also possible to investigate transit times by describing kincmatical paths in the potential 
region [12-14] or introducing a new degree of freedom, like the Larmor clock [9] and time-oscillating 
barriers [15]. Recently, a time operator has been considered [16-18], which is canonically conjugated 
to the energy operator. From this time operator it is obtained a tunneling time that is equal to the 
average over momentum components of the phase time. 

Although the majority of the discussions about tunneling times have been based on the non- 
relativistic Schrodinger equation, some recent works have attended to particular features of relativistic 
tunneling particles, such as superluminality [19-22] and existence of the Hartman effect for relativistic 
potentials [23-26]. The solutions of the Dirac equation are oscillatory in the diffusion, E > V + m 
[27,28], and Klein, E < V — m [29], energy zones, where V is the barrier height and m is the particle 
mass. The phenomenon of multiple peaks occurs in the diffusion zone [30]. The Klein region, in 
its turn, involves pair production and dynamical localized states [31]. The tunneling energy zone, 
max{m, V — m} < E < V + m, is characterized by evanescent solutions [26] . In this paper, we aim to 
present an analytic and numerical study of the phase time for one-dimensional relativistic tunneling 
in the opaque barrier limit. We choose an incoming spectrum and vary the barrier height in such a 
way that the momentum distribution remains at evanescent zone with above potential energies, i.e. 

max{m, V } < E < V + m , (1) 

whence we observe that energy components with E < V a correspond to below potential states whose 
interpretation is still subject to discussions [26,27,31]. 
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In section II, we review some of the standard calculations on Dirac tunneling and find an approxi- 
mation for the transmission coefficient in the opaque limit. This approximation is then used in section 
III to obtain a closed formula for the phase time in relativistic tunneling. Such an analytic expression 
corrects the well-known formula obtained by the SPM. The tunneling time is proportional to the 
barrier width for an incoming momentum distribution whose upper limit is very close to the barrier 
height. This result clearly contradicts the Hartman effect. Our analytical expression and the SPM 
formula coincide for higher potentials. To support our analytical study, we also present numerical 
calculations. Our conclusions arc drawn in section IV. 



II. OPAQUE BARRIERS AND FILTER EFFECT 

Consider a relativistic spin one-half particle of mass m moving along the z-axis in presence of a one- 
dimensional electrostatic potential whose height is V in the region < z < L and zero elsewhere. 
The particle dynamics is described by [h = c = 1] 

id t y(z,t)= lo [m-i l3 d z ]V(z,t) + V(z)y(z,t) , (2) 

where 7 and "f 3 are the Pauli-Dirac matrices defined by 

/ 1 \ , ( 0-3 

7o = ( -1 ) and 73 = { -a 3 



Let us consider as incident wave packet 



with 



% no (z,t) = N f " dp g(p)u(p) e^^-Et) ^ 



g(p)=exp[-(p-p o ya 2 /A} and u(p) = ( 1 — 

a + m 



(3) 



N is & normalization constant and a is related to the spatial localization of the incoming particle. In 
this paper, we choose 



max | , yV~~m 2 j <p m and p M < \JV +2mV 

to guarantee that the energy spectrum remains in the tunneling zone, Eq. (1). This avoids Klein 
and diffusion phenomena, as well the presence of below potential evanescent solutions. By solving 
the Dirac equation (2) and imposing the continuity of the wave function at and L, we find for the 
transmitted wave packet [26-28] , 

% tA (z,t) = N P dp g{jp)T{jp,L)u{p) e^ z - E ^ , (4) 

■'Pm 

with 



T(p,L) 



= e~ lpL 



cosh(gL) — i — — s'mh(qL) 

qp 



(5) 



where q = y/ m — (E — V ) 2 . To summarize our notation, wc shall denote by g m M and £' m M the q 
and E functions respectively calculated in p m and p M - 

For thin barriers, the transmitted momentum distribution is very similar to the incoming distri- 
bution. Consequently, the transmitted momentum is still centered in p„. For very thicker barriers, 
the mean value of the transmitted momentum, (p) T , tends to p m , see Fig. 1. By increasing the barrier 
width, the potential acts as a momentum filter. Observe that the filter effect is more evident for 
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momentum distributions whose upper limit coincides with the upper limit of the tunneling zone, see 
Fig. lb. In the opaque limit, mL 3> 1, we can approximate the expression for T given in Eq.(5) by 



T(p,L) 



mV 



qp e- qL e - lpL+lv 



(6) 



where ip = arctan [(p 2 — EV )/qp\. 

Before beginning our discussion on the phase time formula obtained by using the stationary phase 
method, let us calculate, in the opaque limit, the average momentum 



(p)t = 



I'm 



dpp \\g(p)T(p,L)u(p)\\' 



dp \\g(p)T(p,L)u(p)\\ 



By using the approximation (6), we find 



(p)t 



dpp 3 q 2 g 2 (p)\\u{p)\\ e 



-2qL 



dpp 2 q 2 g 2 (p)\\ u(p)\\ e 



-2qL 



Observing that p dp = E q dq / (V — E) 1 we can change the variable of integration from p to q, 



(p)t 



dqp 2 q 3 [E/(V -E)}g 2 {p)\\u{p)\\- e 
Hp) 



-2qL 



dppq 3 f(p)e 



-2qL 



(7) 



(8) 



(9) 



The filter effect suggests expanding the factors which appear in the integrands around q M , 

{p 2 q 3 )m + (p 2 q 3 ) q ,M (q - 9m) + (q- q u f 



P 2 q 3 



pq 3 = {pq 3 ) M + (pq 3 ) g ,M (q - Qm) + O (q-q M y 
f(p) = fu + fq,u (q ~ 9m) + O (q-q M ) 
After algebraic manipulations, we find 



1 + 



(p 2 q 3 ) q ,M -Pm (pq 3 ) q , 



(p 2 q 3 )m 

Pq.M ipq 3 )m +Pm {pq 3 )q,u ~Pm (pq 3 )q, 



(p 2 q 3 )i 



p* 



2p M (E M -V B )L 



(10) 



Thus, for opaque barriers, the momentum distribution is sharply peaked in the neighborhood of p M . 

The standard phase time formula is obtained by calculating the space-time points in which the 
phase tp is stationary. The maximum of the wave packet is found by imposing that the derivative of 
the phase calculated in (p) T is equal to zero. For opaque barriers, 



{pz-Et- P L + ip) PiU =0 . 
Consequently at the edge of the barrier, z = L, we have E pu r SPM 



Finally, 



Pu<k 



(11) 



We observe that, for E M — > V a + to, t spm — > oo. In the next section, we shall overcome this ambiguity 
by proposing a different method to calculate the phase time. It is based on the analytical calculation 
of the probability density by mean of the approximation (6). We then compare our phase time formula 
with Eq.(ll) and discuss in details the validity of the Hartman effect. 



3 



III. REVISING THE PHASE TIME FORMULA 



In this section, we propose an analytic method to obtaining the phase time. This method is essentially 
based on the search of the time for which the electronic density is greater at the far edge of the barrier, 
z = L. This means taking the derivative of the electronic density with respect to time and finding 
when it is equal to zero, 







(12) 



The time solution of the previous equation will be indicates by r. Using this time, it is possible to 
introduce a tunneling (or transit) velocity, defined by v tm = L/t. 

Numerical solutions, t num , of Eq.(12) were calculated to an incoming gaussian distribution, g(p), 
characterized by 

ma = 10 and (p n ,,p ,p M ) — (0, |, 1) m . 

For a potential barrier with height V = E M — m(w .414 m), we find that the tunneling velocity tends 
to a constant value when we increase the barrier width L, see Table 1 and Fig. 2. We can also see that 
for increasing potentials, for example V = 0.46 m(> E M — m), the transit time tends to a constant 
value. This time is in agreement with the phase time obtained from the stationary phase method, 
Eq. (11), see the last row in Table 1. 

Let us now begin our analytical discussion of Eq. (12). For opaque barriers, we can use the appro- 
ximation given in Eq. (6) and, due to the filter effect, develop the phase around k M (or q M if we intend 
to change the variable of integration from k to q). Using the expansions up to the second order, 



(p = ip M + ip q>M (q -q M ) + V qq , M (q - <7 M ) / 2 + (l-Qu) 

2 r : 

E = E M + Eq >M (q - q M ) + E qqM (q - q M ) I 2 + O y (q - q M ) 
and changing the variable of integration into p = q — q M , we obtain 



Vtra(L,t) 



2 %. £ M«M 

mV (E M -V ) 



- iE M*-«M L 



with 



S(t) := 



dp (p 



e pL exp 



(<P 9 ,m - E qM t) p + 



S(t), 



(Pqq,M 



E, 



(13) 



(14) 



(15) 



At this point, it is convenient to introduce the quantities a and j3 defined by 

V - 2E M q, 



a = a t - a 2 t = ip q<M - E q ^ M t 



Pi P°2 ^ — ^qq.M E qq M t 



p M (E M -V ) V ~E M ' 

= g M [ m 2 K + E M (2E M V - 2 E* M -V*)] 
2p 3 M (E u -V ) 3 



Hv -E u y 



t (16) 



The observation that for q M — > the main time contribution comes from the beta term and that for 
increasing time this term (which is coupled to p 2 ) becomes comparable to the a term (proportional 
to p), suggests to consider the following approximation for the exponential which appears in Eq.(15), 

exp{ i [ ap + (3p 2 } } ~ 1+i (f3p 2 + ap) - \a 2 p 2 ~ \ /3 2 p i - a/3p 3 . 

This means that in the calculation of the integral Eq.(15), we shall find integrals of the form 



s(n) := 
Finally, 

\s(t)\ 2 



dp ( P ^ 2 + 2q M p» +1 +q 2 M p") e~" L 



(n + 2)! 



+ 2 (n + 2)(n+l) 



(17) 



1 2 2 

■ a p 



^^p'-a^p 3 +%{pp 2 +ap)\ 



| s(0) 2 ^ y 2 , 

s 2 (0) + a 2 \s 2 {\) - s(0)s(2)] + 2 a(3 [s(l)s(2) - s(0)a(3)] + /3 2 [s 2 (2) - a(0)«(4)] . (18) 
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Our phase time formula is obtained by taking the derivative with respect to time and imposing that 
it is equals to zero, i.e. 

a a t [s 2 (l)-s(0)s(2)] + (aft + a t ft [s(l)s(2) - a(0)s(3)] + ft [s 2 (2) - a(0)s(4)] = . 

Observing that ftft [s 2 (2) — s(0)s(4)] <C a^a^ [s 2 (l) — s (0)s(2)], we find the following analytical ex- 
pression for the transit time, 



= a^i [s 2 (l) - a(0)a(2)] + (aift + fta 2 ) [s(l)s(2) - s(0)s(3)] 

Tana [ S 2 (l)- S (0)5(2)] + 2a 2 /3 2 [s(l)s(2) - a(0)s(3)] + ft 2 [s 2 (2) - S (0)s(4)] 

In the limit </ M — > (which implies a 2 — > 0), the previous expression simplifies into 



(19) 



Ms(l) a (2)- S (0)s(3)] _ 2(2m + Fo) (3! 4! - 2! 5!)/r _2 /~~2m r , 9m 
ANA ft [ S 2 (2)- S (0) S (4)] VK 2 +2mK (4!4!-2!6!)/L 10 9 V V„ (M) 

The SPM result is recovered for q M values for which the time dependence in a is not more negligible. 
In such a limit, the /3 term plays no role in the calculation and 

a ± _ 2 E M - V a 

Tana r — . l^-U 

a 2 Pm <7m 

The phase time formula (19), the formula obtained by the stationary phase method (11) and the 
numerical simulations are shown in Fig. 3. Our phase time formula is in excellent agreement with the 
numerical analysis. 

Before concluding this section, let us make some observations. In the case, q M — 0, the spatial 
phase exp[i k(x — L)] can be approximated by using 



V + m 



VV (V + m) ° q 2 + O 

2my/V (V +m) 



This implies that the spatial dependence has to be included in the (3 term (which contains the terms 
in q 2 ). To reduce the maximum of S(t), sec Eq. (18), of a factor 1/e, observing that the time (20) 
guarantees a 2 ~ aj3 ~ /3 2 , we have to consider x — L ~ r ANA . Consequently, the spreading of the wave 
packet in configuration space is proportional to the barrier width. 

It is also interesting to note that in the non-relativistic limit, i.e. E -C m and V a -C m, the transit 
velocity L/t ana tends to 4.5 \JV a j2 m, which is clearly subluminal. 



IV. CONCLUSIONS 



In this paper, we have investigated the tunneling through opaque barriers for relativistic particles. By 
an analytical study, we have found a closed formula for the phase time in agreement with numerical 
simulations. In order to avoid "negative energies" in the potential region, the incoming momentum 
spectrum has been restricted to evanescent zone with E > V . We have assumed that tunneling 
time r is defined as the instant of maximum probability of finding the transmitted particle at the 
barrier edge, that is to say, the peak of the wave packet along the time for fixed z = L. Due to 
wave packet spreading, this measurement differs from mapping of the peak dynamics along the z-axis. 
Nevertheless, numerical calculations, which can be readily performed, guarantee that such variation 
is not relevant. An approximation on the transmission coefficient has allowed us to obtain a closed 
expression for r. This formula has exhibited excellent agreement with the numerical calculations and 
has shown cases in which the stationary phase method gives a satisfactory approximation. The most 
important prediction of our analytical formula concerns the validity of the Hartman effect. It does 
not hold for incoming distributions whose maximum energy value coincides with the upper limit of 
the evanescent zone. In this case, the phase time increases linearly as function of the barrier width. 
In this limit the result obtained by the SPM becomes meaningless since its formula gives an infinity. 
For incident wave packets with momentum integration truncated before V + m, we have found a 
tunneling time which is independent of L in the opaque limit, as predicted by the Hartman effect. 
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Supcrluminal transit velocities appear for small barriers, see Fig. 3. This phenomenon surely 
deserves further investigations. Nevertheless, it has to be underlined that there is no causal connection 
between the peak of the incoming wave packet and the peak of the transmitted wave packet. One 
can consider energy components and barrier heights such that new tunneling features are introduced, 
namely negative E — V values and a Klein zone that borders the evanescent one. 

The numerical simulations presented in these paper are based on the numerical calculation of the 
gaussian convolution of the barrier stationary solutions. Braun, Su and Grobe [32] have proposed 
a numerical approach, based on the split-operator technique, to solve the time-dependent three- 
dimensional Dirac equation. In view of the results presented in this paper, it could be interesting to 
revise tunneling phenomena by using the BSG approach. 
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v o = E m ~ m 


V = 0.46 m 


771 L 




^num/^ana 




Tnum/ Tana 


50 


23.62 


0.8807 


6.520 


0.9448 


100 


51.25 


0.9553 


7.339 


0.9906 


150 


77.87 


0.9677 


7.552 


0.9963 


200 


104.3 


0.9719 


7.650 


0.9981 


250 


130.6 


0.9739 


7.706 


0.9988 


300 


156.9 


0.9749 


7.743 


0.9992 


350 


183.2 


0.9755 


7.769 


0.9994 


400 


209.4 


0.9759 


7.788 


0.9996 


450 


235.7 


0.9762 


7.803 


0.9997 


500 


261.9 


0.9764 


7.815 


0.9997 


mr SPM 


oo 


7.918 



Table 1: Tunneling times are listed for an incoming gaussian distribution with ma — 10, p m = 0, 
p a = m/2 and p M = m varying the barrier width L. The barrier height plays a fundamental role in 
their characterization. When V + m is greater than E M , for mL 1 the tunneling time tends to a 
constant value. In this case, the standard formula obtained by the SPM, Eq. (11), represents a good 
approximation. On the other hand, if H- m = E M , i.e. q M = 0, the tunneling time is proportional 
to L. The agreement between analytic, Eq. (19), and numerical data is impressive. 
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Figure 1 




0.2 0.4 0.6 0.8 p/m 



Figure 1: Transmitted momentum distributions, G T (p) = g(p)T(p, L)u(p), are plotted for different 
barrier widths. Continuous line represents the incoming gaussian spectrum with p m = 0, p„ = m/2, 
p M = m and ma — 10. The barrier heights are chosen equal to V = E M — m 0.414 m in (b) and 
V„ = 0.46m in (a). Filter effect is evident in both cases. 



Figure 2 




5 10 15 20 25 30 35 40 45 mL 

Figure 2: Tunneling velocities are plotted as function of mL. They have been numerically obtained by 
the expression (12), with p m = 0, p = to/2, p M = m and ma = 10. Actually, iiV +m is greater than 
E M , then w tun increases linearly with the barrier width (Hartman effect). However, if the incoming 
energy distributions reaches the upper tunneling zone, E M = V + m, the velocity tends to a constant 
value. Numerical tunneling times for L > 50/m are given in Table 1. 
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Figure 3 



mL = 100 
ma = 10 



E 



V2 



m 



SPM 

ANALYTICAL 
NUMERICAL 



0.42 



0.43 



0.44 



0.45 



Vo/m 



Figure 3: Analytical (ANA), numerical (NUM) and SPM tunneling times arc plotted as function 
of Vo/m. The incoming momentum spectrum is characterized by p m = 0, p = to/2, p M = to, 
and ma = 10. The agreement between r ANA and r NUM is excellent. The SPM represents a good 
approximation for higher barriers. 
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